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In this work, we consider a diffusive two-species d— dimensional model and study it in great 
details. Two types of particles, with hard-core, diffuse symmetrically and cross each other. For 
arbitrary dimensions, we obtain the exact density, the instantaneous, as well as non-instantaneous, 
two-point correlation functions for various initial conditions. We study the impact of correlations in 
the initial state on the dynamics. Finally, we map the one-dimensional version of the model under 
consideration onto a growth model of RSOS type with three states and solve its dynamics. 
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I. INTRODUCTION 

Stochastic reaction-diffusion models have attracted recently much interest in (see e.g. Jl],§] and references therein). 
In the last decade, the latter appeared (directly or via mapping) as models for traffic flow H, kinetic biopolymerization 
j|, reptation of DNA in gels H, interface growth ||J^], etc. 

In this context, simple symmetric (SEP) ^ and asymmetric exclusion processes, in one dimension, (ASEP) jl],|]|]] 
play a particular role because of their relationships with integrable quantum spin systems (Heisenberg chains) and 
because of their connection with the KPZ equation |icj |, directed polymers in random media and shock formation 
(see e.g. jl| and references therein). These models have been extensively studied and the ASEP with open boundary 
conditions, as a simple driven diffusion model, exhibits a rich dynamical behaviour involving different nonequilibrium 
phase transitions in the steady states. They can be studied exactly on the basis of the so called Matrix Approach 
(MA), an algebraic approach based on an ansatz for the probability distribution which is related to the integrability 
of some quantum spin chains. This approach provides the full solution of the ASEP (and also the SEP) model, 
including the full phase diagram, density profile and, in principle, any equal time correlation functions. Though, only 
few explicit results are known jHJ about the dynamical correlation functions, much work has been done on the static 
properties. 

The MA has been generalized to solve the stationary states of one-dimensional models with several species jl2| and, 
recently, a first-order phase-transition in some models has been found |l3| (see also jl4| where different results were 
obtained, independently, for the same model). 

The lack of exact results for the dynamics of multispecies models |l5| , |l6| (in particular in dimension d > 1 see also 
and references therein), has motivated us to study in some details the dynamics of a two-species model, which 



is related to the models introduced by Arndt et al. in [h_2 13 1 . We compute explicitly, in arbitrary dimensions, the 



density, and the two-point instantaneous/non-instantaneous correlation functions. We then exploit the exact results 
for the two-point correlation function to study a RSOS-type J6|J7| JT^JT9|] growth model. 

The paper is organized as follows: In section II, we present the general stochastic formalism within which we will 
work. In section III, for the model under consideration, we compute in arbitrary dimensions the density for various 
initial states and in presence/ absence of initial correlations. In section IV, for a translationally invariant version of 
our model, we evaluate the instantaneous two-point correlation functions in arbitrary dimensions. In particular, in 
one-spatial dimension, we assume both cases where initial correlations are absent /present. In section V, we introduce 
and solve a growth model of RSOS-type "with three states" . This analysis is carried out for different initial states 
(correlated and uncorrelated) . Finally, in section VI, we calculate for systems with uncorrelated (but random) as well 
as correlated initial states the non- instantaneous two-point correlation functions. 
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II. THE FORMALISM AND THE MODEL 



Consider an hypercubic lattice of dimension d with TV sites (N = L d ), where L represents the linear dimension 
of the hypercube, and periodic boundary conditions are imposed. Further assume that local bimolecular reactions 
between species A and B take place. Each site is either empty (denoted by the symbol 0) or occupied at most by 
one particle of type A (respectively B) denoted in the following by the index 1 (respectively 2). The dynamics is 



parametrized by the transition rates where a, (3, 7, S = 0,1, 2: V(a, (3) 7^ (7, 5) , T 
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For example the rate corresponds to the process BB — > AB, while conservation of probability leads to 
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The state of the system is determined by specifying the probability for the occurence of configuration {n} at time 
t. It is represented by the ket \P{t)) = J2{ n } P({ n }^)\ n )' where the sum runs over the 3 N configurations (N — L d ). 
At site i the local state is denoted by the ket \rii) = (10 0) T if the site i is empty, \rii) = (0 1 0) T if the site i is 
occupied by a particle of type A (1) and \rii) = (0 1) T otherwise. It is by now well established that the Master 
equation governing the dynamics of the systems can be rewritten as an imaginary-time Schrodinger equation: 



d_ 

Of 



P(t))=-H\P(t)}, 



(1) 



where H denotes the Markov generator, also called stochastic Hamiltonian, and is in general neither hermitian nor 
normal. Its explicit form is given below. We also introduce the left vacuum (x\ which is defined by 



(x\ 



£<Wl 
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Probability conservation yields the local equation (stochasticity of H) 

(Xl#m,m+e° = 0, 

where e a denotes the unit vector in the direction a (1 < a < d) and m designates a point of the hyperlattice 
labelled with help of its d components: m = (mi, . . . , rrid) ■ 

In this work, we assume that there are only symmetric nearest-neighbour jump-processes. A particle A (respectively 
B) can jump, with rate TjJ = T^J > (respectively = T^q > 0) to an adjacent site (in the d directions) if the 
latter was previously empty. Such processes are symbolized by the "reaction" A0 < — > 0A (respectively, B0 < — > 0£? 
). In addition we assume that when two different particles A and B are adjacent, they can cross each other with rate 
T\\ = r^j > 0. This processes are schematized by the reaction AB < — > BA. 

The local Markov generator corresponding to this system and which acts on two adjacent sites m and m + e Q reads 
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(3) 



where the same notations as in reference [|21|,|17| have been used. Probability conservation implies that each column 
in the above representation sums up to zero. 

Locally, the left vacuum (x\ has the representation (%j = (111)® (111) 

The action of any operator on the left-vacuum has a simple summation interpretation. This observation will be 
crucial in the following. Below we shall assume an initial state \P (0)) and investigate the expectation value of an 
operator O (observables such as density etc.) : (0)(t) = (x\Oe~ Ht \P(0)) 

From (|^) , we can compute the equations of motion of the density and of the two-point correlation functions [ fl7| . 
For the density, we have: 
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For the derivation of the equation of motion of two-point correlation functions we would proceed similarly, as in 
jL7|. One should however pay attention to distinguish the case of the correlation function of adjacent sites from the 
general case. 

In general, when T^J, rJJ, T^q, r^, T^, I^f are independent parameters the equations of motions of the multi- 
point correlation functions constitute an open hierarchy and the dynamics is not soluble. The stationary states of 
such systems have been studied in fll2]| by Arndt et al. Recently it has been shown [fl3f with help of quadratic algebra 
techniques [ ^3|Jl2| and numerical means that an asymmetric version of this model exhibits a first-order phase transi- 
tion, in its stationary state, when T^J = Tq° = rS^ = = = 1 and Tjl = q. The steady-state of the density of 
the same model has also been studied independently by Rajewsky et al. who obtained different results: in Jli) , 
the authors, argued that there is no phase transition from the "mixed phase" to a "disordered phase" . 

Here we assume 

pOl _ plO _ r 02 _ F 20 _ r 21 _ pl2 _ p /r\ 
1 10 — 1 01 — 1 20 — 1 02 — 1 12 — 1 21 — 1 ! \°) 

which guarantees that the equations of motion of the correlation functions close in arbitrary dimensions. 

From now on we focus on the soluble model described by (^|) with equal rates, according to the solubility constraints 

®- 

Before studying statistical and dynamical properties of this model, let us comment on its solvability. In the 
single species reaction-diffusion models, the solvability inherent to the closure of the hierarchy was explained in the 
framework of the duality transformations. In fact, it has been shown that the spectrum of the single-species stochastic 
Hamiltonian (with the solubility constraints) is identical to the spectrum of an anisotropic spin-1/2 Heisenberg 



quantum Hamiltonian Hxxz in a- magnetic field |20|. As shown in |2l| , the situation is quite different for the 
multispecies problem and a general, comprehensive and unified understanding of the formal solubility is still lacking. 
However, for the model under consideration here, it has been shown |22[ that the stochastic Hamiltonian (|J), can be 
mapped, via a similarity transformation, to an exactly integrable quantum spin-1 model introduced by Sutherland 



III. EXACT STUDY OF THE DENSITY 



In this section, we study the density of the system, in particular, when translation invariance is broken (the initial 
density is non-uniform) and when correlations in the initial state are present. 

It follows from (|J) that the density of species j € (A, B) at site m, labelled with its d components (m = 
(mi, . . . , wid)), obeys to the following linear differential-difference equation ]T^ ]: 

j t {nL)(t) = -2rd{nL)(t)+rJ2({nL + e«) + (<-«*)) (6) 

a— 1 

We first consider the situation where particles are initially non-uniformally distributed. Namely, we assume that 
particles of type B are located in the region of space L/2 < x% < L, . . . , L/2 < Xd < L while particles of type A are 
initially confined in the region < x\ < L/2, . . . , < Xd < L/2 (we assume that L is even). Within each of the two 
regions, particles of each type are distributed uniformally with respective densities p_b(0) and /^(O). Solving (^J) for 
this initial condition, we find: 



(t)) = p A (0)e- 2rfr *^«,(0)) J] I ma -m' a (2Tt) = PA {Q)e- 2drt Y. lit 9 ( ? - m « ) W^(2rt)] 



. - 2 

m a—1 



d 



= Pa(0)Y[\ e- 2rt I m ^ m/ j2Tt) , (7) 

a=l \0<m,' a <% j 

where I n (z) denotes the usual modified Bessel function. We have also introduced the Heaviside function Q(x) = if 
x < and 8(x) = 1 if x > 1. 

n 2 n 2 

With help of the asymptotic expansion of the Bessel functions and approximating ^ n e~~ « J c?n Q e _ «r, we 
obtain from <M) the long-time behavior of the density at site m: 
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where erf(z) denotes the usual error function. 
Similarly, we have for the density of particles B: 

d I \ d 

<»£(*)} = PB (0) J] E e- 2rt I ma . m ,j2rt) ] « p B (0) [] 

Q=1 \T <m a^ L 
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We now pass to the case where the distribution of particles for each species j E (A, B) is given by (Kj denotes a 
real dimensionless constant) an initially correlated distribution: 



(<)(0)=ft 



(o) ( n <w° + *j n 
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(10) 



The exact densities then read 111'; 



(<>(*) = Pj (0) 



n (e^/ m ,(2R))+ K , n 
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For 7^ 0, in the limit m ~ L 3> 1 and ^> 1, with cr = m/L and u = When cr = 0(1), then 



we obtain 



(11) 



cr L m 

4rt — 4Tt ' 



n — ' 
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(l-7 i )(4-uer 2 rt)Ti/ 2 
n d (l+e— 2 ")C(7») 

lli=i ( 7r rt) 1 /2 

(e~ g2 " ln(4MoTt)) d 
(7rrt) d / 2 



if < 7i < 1 
if 7ii > 1 
if7i = 1, 



(12) 



where = X)fc>i ,v > 1 is the Riemann zeta function. 

It follows from these results that for 7$ 7^ 1, the density decays as a power-law of time. However, notice that when 
initial correlations are strong (i.e. < 7^ < 1), the algebraic decay of the density is non-universal (it depends of 7$). 
When initial correlations are weak (i.e. Tj > 1), the algebraic decay of the density is universal. Hence the case where 
7t = 1 is marginal and there are logarithmic corrections to the universal behaviour. 

For initial states decribed by (|l0|), with 7^ = and Kj 7^ 0, then we have 



-da u 



r(i/2)-r(i/2,«) 
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(13) 



Therefore the dimensionality has a non-trivial effect: when d < 2, the densities decays as f d / 2 . Otherwise, when 
d>2,{ni a )(t)~t-\ 

On the other hand, when Kj = 0, the initial density of species j vanishes on the hypercube except at the origin, 
where its value is Pj(0). In this case, the limit considered above, yields 



(<)(0) 



~do 2 u 



{Ai:Tt) d l 2 



(14) 



Notice that because of conservation of the number of particles, in the translationally invariant situation, we simply 
have 



PA(t) = pA(t — 0) = PA ;PB(t) = pB(t — 0) = pB 



(15) 
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IV. INSTANTANEOUS TWO-POINT CORRELATION FUNCTIONS FOR TRANSLATIONALLY 

INVARIANT SYSTEMS 



In this section we compute exactly the two-point correlation function for translationally invariant systems, in 
arbitrary dimensions for different initial states. 

The equations of motion for the connected correlation functions C l J (t) = C!f r (t) = {n\n{^} (t) — piPj — (n-o n m-i) (t) — 



PiPji ihj) G (A, B), read, with the notation: r = (ri, . 

n d 
a 



, rd) = m — 1. where a = 1, . 



-C?{t) = -ATdCi\t) + 2r£ a«(t)+C^W , I Ml > 2 



dt 

d_ 

dt 

d_ 
dt 



C%(t) = 2T 
#'(«)= 



<&.(*) + £ {<£-««'(*) + (*)} - (2d - W-W 



(16) 



Solving the latter, we have, using known properties of modified Bessel functions (see appendix, (|5§|)): 

d d 

C?(t) = C l J(0)e- Mrt l[ I ra ^(ATt)+C^(0)e- idTt [] 7 r „(4rt) 

a— 1 a— 1 

- f* dre~ AdTT C% (t-r)(J^- 4dvj J] 7 ra (4Tr) - MT J* dTe~ idrT C% (t - r) JJ I Ta (4Tr) 

+ 2r / rfTe- 4drT ECe«(t-T)[7 r£S+ i(4rT) + / ra _i(4rr)] J] / rQ ,(4rV) (17) 

•'° a=l a'=jta 

Restricting the solution to one-spatial dimension, with r = \m — l\ > 0, the Laplace transform yields (i,j£ (A, B)) 

C\\t) = e- 4rt Y, C r> (0) {^(4T*) + Ir'-i(4rt)} (18) 



r'>l 



and more generally, 



Let us now consider one-spatial dimension and assume that the initial correlations are given by 

C l r {0) = Kir'" 1 , vi > 0, I e (AA, BB, AB) 



(19) 



(20) 



We discuss the case |k/| > while the case m = corresponds either to the situation where no particle is present 
on the lattice initially, or, when all sites of the lattice are occupied by particles of species i (or j). An alternative 
is that the system would be initially in its steady-state. When a single species is present initially, say species A, we 
recover the known problem of symmetric diffusion of hard particles A + <-> + A. When the lattice is full (or empty) 
initially, no dynamics takes place. A single-species one-dimensional process A + A <-> + with a correlated initial 
state as in ( pGf ) has been studied in |Q . 

Again we can infer the asymptotic behavior of the two-point connected correlation functions in the limit Tt 3> 1 

With U EE £ < 00 0. 

It is useful for the sequel to introduce the definitions of the auxiliary functions |17| : 

_ (r(i^) + r(i=* a*v) - r(i=* «(i - of) - r(±^, «(i + a) 2 )) 



Fi{v,a,vi) 
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We distinguish two regimes 



i) For r < L, with r^/STt < 1 and a = f 



L > 



\/4¥(8rt)"i/ 2 

Kl(2C(w,)+(8«rtr|i)" 



(8wrt) 1 /2 

Kj In (8r-»(l-o- Q )t) 
(87rrt)!/2 



if < n < 1 

if n > 1 
if zy ; = 1, 



(23) 



ii) For r» l,r = ah ~ L, we have 



(i+e— 2 ")c(^)+((i-^)/^)(8«r ( r 2 t)- : 



(87rrt) 1 /2 



Kje - * 7 11 ln(81><Tt) 



if < vi < 1 

if l/jr > 1 

if = 1, 



(24) 



As for the density, it follows from these results that for vi ^ 1, the (connected-)correlation functions decay as a 
power-law of time. When initial correlations are strong (i.e. < U\ < 1) the power-law decay of the correlation 
functions is non-universal . In contrast, when initial correlations are weak (i.e. vi > 1), the algebraic decay of the 
correlation functions is universal. The case where v\ = 1 is marginal and logarithmic corrections to the universal 
behaviour arise. 

In arbitrary dimension (d > 1), we consider a translationally invariant random but uncorrelated intial state, 
described by: 



Cr(0) =Ki,le {AA, BB, AB), 



(25) 



where, as above, m ^ 0. 

The asymptotic behaviour (L,Tt^> 1 with v = L 2 /8Tt < oo and a a = r a /L) of the connected correlation functions 
is then @: 



C' T (t) 



Kl 



I 1 + T£rt) d Tit=i -T 7 !,^, cr Q , = 0) si r a < L 



ki (l+ift) n Q =i- ?r 2, Q («,c r a,^ = 0) sir a ~£»l, 



(26) 



where the quantities T\. a and Tn, a are obtained, respectively, from ( ^l|) and ( p2|) on substitution of cr by a a = r a /L. 



V. THE MAPPING ONTO A RSOS GROWTH MODEL WITH "THREE STATES" 

In this section, we introduce a growth model of RSOS type (Random Solid on Solid ) with "three states", by 
exploiting a mapping of the one dimensional model studied in the previous sections. 

Let us briefly recall that the RSOS growth models are e.g. useful to describe the spatial fluctuations of the (one- 
dimensional, of length L) interface location in the magnetization profile between coexistent phases in two-dimensional 
models of ferromagnets, such as in the zero field planar Ising model In such models, at zero temperature, every 
path minimizing the energy of the systen is a sequence of L binary numbers Uj = 0, 1 with j = 1, . . . ,L. The stochastic 
variable Uj has the value nj = if the j th segment of the interface steps upwards (in an angle of 7r/4). The value 
iij = 1 corresponds to the case where the segments steps downwards with an angle of 7r/4. The quantities nj = 0, 1 
can be interpreted as occupation numbers relating the interface height hj according to hj — /ij-i = 1 — 2rij In this 
case the displacement Sh r (t), at time t, of the segment of the interface from the sites fa to fa > ji, with r = fa —fa > 
is given by Sh r (t) = 5^I =1 (1 — 2nk(t)) and thus \Sh r+ i — Sh r \ = 1. 

Here we consider an extension of the above model. We consider that the configurations minimizing the energy (at 
zero temperature) are of the form {xAn± 4- xsnf , xa^ + xsnf , . . . , £a^l + xsnf }, where the discrete stochastic 
variable XAnf + XBTif can take three values. The case XAnf + xsn^ = again corresponds to the situation where 
the j th segment steps upwards with an angle of n/4. The case XAnf + xgnf = xa (respectively, XAnf + XBnf = xb) 
describe the situation where the j th segment forms an angle Arctan(l — xa) (respectively, Arctan(l — Xb)) with the 
horizontal. When xb — (respectively xa — 0) and xa — 2 (respectively xb = 2), we recover the above-mentioned 



G 



two-state RSOS growth model. In this sense the growth model which we study hereafter is a "three-states" extension 
of the usual RSOS model [p|j7| |l^ , [l^] . In addition the mapping with the diffusive model is clear: the presence of 
particle of species A (resp. B) at site j translates in the language of the growth model with the fact that the related 
segment of the interface forms an angle Arctan(l — xa) (respectively, Arctan(\ — Xb)) with the horizontal. In this 
picture, the jumping of the diffusive particles correponds to the fluctuation of the orientation of the related segments 
of the interface. 

We consider a translationally invariant system; the displacement of the (one-dimensional) interface, at time t, from 
the sites k\ to &2 > &i , with r — k% — k\ > is given by hk 2 (t) — ^fei (*) = Sh r (t) , where 

T" 

Sh r (t) = Y, (1 - XAni - x B n%) (27) 

m— 1 

Therefore, in the model considered here, we have max\5h r+ i(t) — Sh r (t)\ = max(\x A — 1|, \xb — l|j \xa + x b — 
instead of the usual constraint \8h r +\(t) — 5h r (t)\ = 1 of the conventional RSOS models. The mean-displacement of 
the interface reads (Sh r (t)) = r(l — XApA — xbPb) thus , if one wants to impose a zero mean-displacement of the 
interface , we have to require that xaPa + xbPb = 1- 

In this section we are interested in the computation of the fluctuations of Shj (t) : 

r r— r'+l 

((Sh r (t)) 2 ) = w 2 (r,t) = ((xa+x b ) 2 - l) r 2 + £ £ [x 2 A (nfnf +r ,)(t) + x 2 B (nf nf +r ,)(t) + 2x A x B (nfn? +r ,)(t)] , (28) 

r'=l 1=1 

where w(r, t) is the physical width of the interface . 

Using the fact that §- t ELlKXn+zXX) = 2F (<«m n m+r+l)(*) - ( n m n m+r)(*)) and El=l E[=l ^("m n m+/>(*) = 
2r ^(rf m nm+r+i)W — ( n m n m+i)(^)) j we obtain the following equation of motion for the width: 

dw 2 (r,t) 



dt 



AY [x\ {C AA (t) - C AA (t)} + x% {C BB {t) - Cr(t)} + 2xaxb {C AB (t) C AB (t)}} , r > 1, (29) 



with w 2 (l,t) = w 2 (l,t = 0) = (xa+xb) 2 + 1 - 2(p A XA+ p B x B )- 
For r > 1, we thus have: 

w 2 (r,t) = w 2 (r,0) 

+ 4r /' dt' [4 {C r AA (i') - C AA (t')} + x\ {C BB (t>) C BB (t')} + 2x a xb {C AB (t>) - C AB (t')}] , , (30) 
Jo 

where 

iu 2 (r,0) = [(x A + x B ) 2 + 1 - 2(p A x A + PBX B )]r 2 

r 

+ 2 £(r- r') [4(^„^ +r ,)(0) + 4<n*n£ +r ,)(0) + 2^^(^n^ l+r ,)(0)] (31) 
r'=l 

With help of the formula @: 4r J * dt' e'^' I n {ATt') = ATte^ 1 (I (4Tt) + h(ATt)) + n ( e - 4r %(4rt) - l) + 

2e~ 4r * Efe=i ( n — k)Ik(ATt), and with the explicit expression of the correlation functions ( |Is| ) and ([H]), we obtain the 
following exact expression for w(r,t): 



w 2 (r,t) - w 2 {r,Q) = 2e" 4rt £ + *|^» + 2x 4: r B C; 4B (0)) 

r'>l 

+r'-2 r+r'-l r'-2 

^ (r + r'-fc-l)4(4rt)+ ^ (r'-fc-r)7 fe (4R)-2^(r'-fc)4(4ri)-7 r ,_ 1 (4n) 



fe=i fe=i fc=i 



(32) 



We will now specifically focus on two kinds of initial states: 

(i) We assume first that the system is initially characterized by an alternating periodic array of particles of type 
A and B. We consider thus an initial state |P(0)) = \ (\ABAB . . .) + \BABA . . .)), with x A + x B — 2 (in this case 
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(6h r (t)) — 0). This initial flat interface leads to pa — Pb — 1/2, and for the connected initial correlation functions, 
we have: C AA (0) = ^ = C BB (Q) = -C AB (0). 

It follows from ( |3l| ) that the initial fluctuations read in this case: w 2 (r, 0) = r (4r + \{xa — xb) 2 — l), for r > 1 
and r even. Therefore, we also have w 2 (l,t) = 3. 

For this initial configuration, the expression ( |l8| ) simplifies and we have C AA (t) = C BB (t) = —C AB (t) = 
— — Io(4Tt). The general expression ([32] ) of the fluctations reads: w 2 (r,t) — w 2 (r,0) — (xa — 
x B ) 2 T J* dt'e~ 4Tt ' 7 (4rt') + ^ r />i(-l) r ' (I r+r '_i(4iy) + J r _ r /(4IY)) . From this expression, using the asymptotic 
behavior of the Bessel functions, it is possible to obtain the long-time behavior (for Tt ^> 1 r 3> 1) of the fluctuations 

w(r, t) 2 - w(r, 0) 2 = (x A - ^) 2 ^(1 + ©((H)" 1 )) 

+ [x A -x B ?Y. f ~4= {e-^-^rf + e -(^')Vs^ } (i + 0((r _ /)- >)) (33) 



r'>l 



Thus for the initial condition considered here, when a; a 7^ sb, it follows from (|33|) that the fluctuations growth as 
w(r,t) « (ri) 1/4 . 

On the other hand, it is known that for the flat interface (or "sawtooth" initial state) , which dynamics is coded in 
a two-state model 5h r = 2m=o(^ ~ % n m)> tne fluctuations growth as ~ (rt) 1 / 4 [^[[T^p. In the situation considered 
here, the fluctuations still grow as w(r,t) ~ (Tt) 1 / 4 and thus the details of the model and its "three-state" character 
only appears through the amplitude (xa — xb ) 2 ■ When xa = xb = 1, the initial configuration corresponds to a 
straight line and, according to (|33"|), there are no fluctuations. Let us also note that when, e.g., xa — 2 and xb = 0, 
the B particles play the role of the vacancies in the two-state model and the model ( p7|) is exactly mapped onto the 
well studied two-state model [j? 18 |] . 



(ii) We now investigate the fluctuations of 5h r (t) in the presence of initial correlations. 
We assume that particles of type A and B are distributed according to Jl8[ 

4(0) = /c;r-" ! , vi > 0, Z e (AA, BB, AB) (34) 

which corresponds, via the mapping (^7|), in the language of the growth model, to an interface with initial fluctuations 
given by ©, where «< +r ,) (0) = C*j(0) + € (A, B) 

Let us define the following quantity: 

f = min(vAA, vbb,vab) (35) 

With help of (p9| ) and (|2(J), (p3|), we can compute th asymptotic expression of the fluctuations for Tt ^> 1 and 
r» 1, with i> = and a = r/i = 0(1) (we assume that x axb 7^ 0), which reads 

SZ v (llrty/* CT ' ^) ~ ^ 2 ( u > CT ' , if < < 1 

u;(r,i) 2 - W (r,0) 2 ~ { V2£k(u)(l - e-°* v ) , if 14 > 1 (36) 

2££ (ln(8rwt) - e"^ 1 ' ln(8IW*)) , if v% = 1 



where the quantities T\ and .F2 have been defined in ( |2l|), ( |22| ) and £(y) is the usual Riemann zeta function. 

We see that in the presence of initial correlations (|34|), the fluctuations are dominated by the smallest initial 
correlation exponent v. Therefore, if v = vaa (resp. v — vbb ), the dominant contribution ( p6|) to the fluctuations 
are the same as for a correlated two-state model where xa = 2 and Xb = (resp. = and xb = 2) where the £? 
(resp. the A) particles play the role of vacancies |Ts| ]. 

From (|3^) we see that initial correlations affect the long-time behavior of the fluctuations of the height displacement 
of the interface: when the correlations are "strong enough" (i.e. < v < 1), thus w(r,t) ~ (Tt)z~ ! i. Conversely, for 
"weak" initial correlations (v > 1) we recover the usual fluctuation exponent: w(r,i) ~ (Tt) 1 / 4 . The intermediate 
case v = 1, corresponds to the marginal behavior where w(r,t) ~ (Li)* vln Ti 

VI. NON-INSTANTANEOUS TWO-POINT CORRELATION FUNCTIONS 

In this section, we compute exactly the non-instantaneous two-point correlation functions for various initial states. 
Similar quantities have already been computed, for some specific single- species models (see e.g. 18 25|- 127|] ) . 
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1 - pa(P) - Pb{0) 

Let us first consider an uncorrelated initial distribution |P(0)) = | Pa(0) 

Pb(0) 



, such that 



(nt(0)nf (0)) = PA {0)5 mA + PA (0) 2 {1 - 8 m ,i) ; (n^(O)nf (0)) = p s (0)<5 m ,i + p s (0) 2 (l - <5 m j); 
(«4(0)nf (0)> = <«m(0K(0)) = p A (0)p B (0)(l - <W) (37) 
We then have : 

(«1W«i 4 (0)>-P? 1 (0) + (pa(0)-p^(0)) n e- 2r *7 roa _ ia (2R) (38) 

a=l...d 



(n±(t)n?(0))=p A (0)p B (0) 



n 



„-2rt ; 



<m Q -Z, 



(2rt) 



(n*(i)nf(0)) 



(39) 



U*)nf(0))=p|(0) + (p B (0)-p|(0)) J] e- 2ri / mo _j2R), 



=l...d 
'2 , 



(40) 



We are interested in the asymptotic behavior (Tt ^> 1 and u = L /4Tt < oo) of the above functions in the two 
regimes: 

i) \m a — la\ = f a ~ L ^> 1, in this case cr Q = r a /L = 0(1). 

ii) \m a — Z Q | e r„ « 1, in this case a a — r a /L = 0(1/L). 

It is worth noting that the autocorrelation functions are obtained in the second regimes (ii). We then have 



<^K(0)> = Pi(0) + + ^(Vt d ) 



(41) 



(«4(*)ni fl (0))=PA(0)p B (0) 



(47rR) rf / 2 



0{\/t d ) 



= «(*)<(0)) 



(42) 



<(*K(0)>=/>S(0) + 



(p B (0)-p|(0))e-E: =1 ^ 



C(l/< d ) 



(43) 



(4^)^/2 

In these regimes, we have a power-law decay of correlation functions € (A, B)), namely, 

(nUt)n{(0)} ~ (rtr^e-^y*- (44) 
Let us now pass to the case where the initial state is correlated according to 
«(0)nf(0)) = Kij Yl (1 5 rafi )\r a \- A Z, r a = \m a - l a \, A§ > 0, (ij) G /Q, > 0,dist{l,m) > 0. (45) 



= l,...,d 



and 



(nj n (0)ni 1 (0))=Pi(0)p J (0)% 
Notice that in one dimension the initial state ( ^o] , ^6|) is translationally invariant and reads: 

(r4(0)nf (0)) = « = | m _ ; |(0)r4(0)) - - 5 r ,o)r" A « + P MkA^ = **5 ^ = ^ 



(46) 



(47) 



This translational invariance which is broken (see (j45|)) in higher (d > 2) dimensions leads to two regimes: 

i) We begin with the one-dimensional case (d = 1), here r = r a = \m — l\. Because of the initial translational 
invariant state, we expect that the non-instantaneous correlation functions only depends on r = m — I, and we obtain 
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(ni(t)nf(0)) = (n?(t)rtf (0)) - ^(0) e 2rt / r (2R) + K AA £ l»T A ^e" 2r '/ r _ r , (2R) (48) 
<»£(*)nf (0)) = ICab J2 \r'\- AAB e-^ Ir _ r , i2 rt) = (ng(t)nf(O)) = (ng(t)nf(0)) (49) 

«(i)nf (0)) = (nf (t)n B (0)) = p s (0)e- 2rt / r (2ri) + /Cbb ]T |r'|- AsB e- 2r 'W(2n) (50) 

Because all the processes in the evolution operator are symmetric (unbiased), there is no drift therefore 
ii) In higher dimensions (d > 2), 
(n±(t)nf(0)) = PA (0)e- 2dFt JJ I m ^ 7< {2Tt) + K AA ]T J] K ~ Q^^^ Im a - m ' a (2Tt) (51) 

a=l...d (m' 1 =jtl 1 ,...,m' d ^tl d ) a=l...d 

(ni(t)n B (0))=lC AB ]T II K-^r A ^e- 2r % ls _w(2rt) = (n^(t)<(0)) (52) 

(m' 1 =£l l ,...,m' d =£l d ) a=l...d 

(n B ,(t)n B (0))= PB (0)e- 2drt ]J I m ^ r< (2Tt) + K BB £ J] \m' a - l a \- A BB e ' 2rt I ma _ m ,j2Tt) (53) 

a=l...d (m' 1 ^Zi,...,mJj^i ci ) a=l...d 

(ng(t)nf(0))=JC BA £ J] |m' ct -/ ct |- A S- e - 2rt / mi _ m ,(2R) (54) 

( ro^T^i,. ..,mj,7^<i) Q=l...d 

We observe that in higher dimensions, because of the initial correlations, the non-instantaneous correlation functions 
no longer depend on \m a — l a \. 

We can express the asymptotic behavior of these non-instantaneous correlation functions in an unified way including 
both d = 1 and d > 2 cases. Assuming that r a = \m a — l a \ <~ \m a \ ~> 1, with r a — <r a L and u = L 2 /4Tt < oo, 
Tt,r 3> 1, the asymptotics reads ((i, j) G (A, B)) 



Kn(*K(0)) 



p t (0)e tt 
(4^) d / 2 U 

v ' a=l... 



+ 0{r 2d ) I , < A* < 1 (55) 



_1-A«V TT 4 U r ( x2t A S / 2_ 

Moreover 

<r4(*K(0)) = (4 ^ )d/2 ^(Oje-S^^. +/ Cij [] « A S) +°( rM )) > A S > 1 (56) 
When Afj = 1, we have the marginal case with logarithmic corrections 



(4 



-L^-^^e-E^^ +ZCy J] ln(4 U a«R)+Or 2d ))^ , A" = 1 (57) 



Again, strong initial correlations lead to an algebraic decay of correlation functions (0 < Ag- < 1), ("mW n i(0)) 
^ , < Af 5 < 1, while for weak initial correlations, we have («mW n i(0)) ~ 

(4 



^ — - — , < Afj < 1, while for weak initial correlations, we have {n t m {t)n\ (0)) <~ ^rtW 2 ' A S > ^' ^ e 

marginal case Ag = 1 is characterized by (n l m {t)n\ (0)) <~ ^"p^ /2 , Ag- = 1 
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VII. SUMMARY AND CONCLUSION 



In this work we studied, by analytical methods, the dynamics of a symmetric two-species reaction-diffusion model 
in arbitrary dimensions. We mapped this model onto an one-dimensional RSOS-type growth model and obtained 
explicit results for the latter. In particular, we were able to compute the density profile for three various initial 
conditions (uniform and non-uniform) in arbitrary dimensions. 

Furthermore, we evaluated, for a translationally invariant system, the instantaneous two-point correlation functions 
in arbitrary dimensions. In one-spatial dimension, we considered the case where initial correlations were present. We 
observed that when the initial correlations are strong enough, they affect the asymptotic dynamics. We also noticed 
a crossover in the dynamics between the case of strong and weak initial correlations. 

We mapped the one-dimensional version of the reaction-diffusion model on a "three-states" RSOS-type growth 
model. Using the exact, instantaneous correlation functions, we computed the exact expression of fluctuations of the 
interface for the latter model. We specifically considered the case of a "flat interface" where the fluctuations growth 
as (Ft) 1 / 4 , as in the corresponding "two-state" growth model: the three-state nature of the model considered only 
appears in the amplitude. We also considered the case where the initial configuration is translationally-invariant, 
random and correlated. We saw that the initial correlations are "strong" , they affect the long-time behavior of the 
fluctuations of the displacement interface. Conversely, "weak" initial correlations, do not affect the dominant term 
and the fluctuations still grow as (Ft) 1 / 4 . This is analog to what happens in the two-state RSOS systems where initial 
correlations affect the long-time behaviour of the width |Q . 

Finally we computed in arbitrary dimensions the exact non-instantaneous two-point correlation functions for initially 
uncorrelated states as well as for cases where correlations were present. Here we again observed the effect of strong 
initial correlations on the dynamics and a crossover between regimes with strong and weak initial correlations takes 
place. 

We conclude this work by addressing an interesting question based on the similarity of the stochastic Hamiltonian 
under consideration with the integrable Sutherland's quantum spin system p^ , p3| . It is known that for one-dimensional 
SEP-model, the relation with the Heisenberg chain has been fruitful to obtain a ( dynamical) matrix formulation of the 
probability distribution, which allowed to solve the (dynamical) density profile for a SEP-model with open boundary 
conditions (particles were injected and ejected from both ends of the chain) ||. The relation of the present model to 
Sutherland's one suggests that a dynamical matrix approach would be possible (in one-spatial dimension) to treat the 
injection and ejection of particles at the boundary (open boundary conditions). 



The support of Swiss National Fonds is gratefully acknowledged. 

IX. APPENDIX: SOLUTION OF THE EQUATIONS OF MOTION OF THE INSTANTANEOUS 

TWO-POINT CORRELATION FUNCTIONS 

Solving the equations of motions ( |l6|) for the instantaneous correlation functions in arbitrary dimensions, for the 
model under consideration in section VI, we obtain the following expression: 
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d d 



<#(*) = C%(0)e- 4drt H I ra (4Tt) + ]T C^(0)e- 4rdt J[ I ra . r , a (4Tt) 





+ 2r f dt'e- 4 ^ 



o 



a^a' — l...d 



H I ra ,(4T(t-t'))(I ra+1 (4T{t 



t'))+Ir, 



x(4T(f 



0)) 
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4dr f dt'e- iTd ^^C%{t')T\ I r {AY{t-t')) 

J a =l 



(58) 



Using the properties of the derivatives of Bessel functions and then integrating by parts, we obtain the more compact 
form (In 
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